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1. Introduction 

Determinations of the bottom quark mass nib have been the subject of a large number of QCD 
precision studies in the past. For a summary we refer to Ref. [1]. The bottom mass is an important 
parameter in numerous theoretical predictions not only within, but also beyond the standard model. 

The data from e + e~ collisions is the common experimental input in many determinations of 
nib, because in particular the region close to the bb threshold and the T resonances of the total cross 
section are very sensitive to the bottom mass parameter. One classic approach is based on the sum 
rule [2] that states the equality of the experimental moment 

P n P = I (1.D 

where R bb - = c(e + e~ — > bb + X)/o pt is the measured inclusive (hadronic) bottom pair produc- 
tion cross section normalized by a pt = 4na 2 /3s, and the corresponding theoretical expression P'J 7 
obtained from an operator product expansion (OPE) in QCD. For not too large n nonperturbative 
power corrections to P'J 1 are suppressed and the theoretical prediction is dominated by the perturba- 
tive QCD (pQCD) result for an external bottom quark pair. Concerning the appropriate theoretical 
formalism bottom mass determinations from the equation F* (nib) = Pn* p differ depending on the 
values for n. We distinguish two classes. 

For n < 3 the theoretical moment P'J 1 is governed by fluctuations at the scale nib- Therefore 
higher order terms in the OPE typically scale like powers of Aqcd / nib and the conventional pQCD 
result is in principle sufficient for a precise determination of the bottom mass. Recent low-« analy- 
ses [3, 4] employ an approximate four-loop pQCD calculation. Ref. [5] uses a related variant of this 
method, where the integration in Eq. (1.1) is only carried out over a finite range and a compensating 
term (according to Cauchy's theorem) is added to the theory prediction. Both Iow-ti approaches 
have the drawback that precise experimental data is currently only available in the region close to 
the production threshold (and for the T resonances) and this deficiency has to be compensated by 
additional theory input in one way or the other 1 , see e.g. Refs. [6, 4] for discussions on the cor- 
responding uncertainties. Power corrections to the additional theory contribution 2 in the moments 
are commonly assumed to be small. 

On the other hand large-« moments, where 4 < n < 10, receive only negligible contributions 
from the energy regions beyond threshold and are dominated by the experimentally well-known 
T-resonances and hence nonrelativistic bound state dynamics. The corresponding sum rules are 
therefore often called nonrelativistic or T sum rules. Due to the nonrelativistic nature of the large-/! 
moments in addition to the hard scale nib the s °ft scale nib / y/n and the ultrasoft scale m^/n emerge 
as relevant short-distance scales and the convergence of the OPE requires the upper limit n < 10. 
The hierarchy between these scales induces sizable terms oe (a s y/n) k , the so-called Coulomb sin- 
gularities, and large logarithms oc (a s ln(«))' in the perturbative loop expansion. The resummation 

'in fact the finite energy sum rule used in Ref. [5] is equivalent to an infinite energy sum rule, if in the range 
above the finite energy limit the theory result is used for the total cross section in Eq. (1.1) and the infinite moment 
integration converges. 

2 related to the fact that the energy integration contour needs to be deformed onto the positive real axis close to the 
finite energy cutoff 
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of the Coulomb singular terms to all orders can be performed within the effective field theory 
NRQCD [7, 8]. Extensions of this framework like the pNRQCD [9, 10] and the vNRQCD [11] 
formalism also allow the systematic resummation of the logarithmic terms. The renormalization 
group improved (RGI) result for the theoretical large-?! moments is expressed as a simultaneous 
expansion in a s and 1 / y/n and schematically takes the form 



k.i 



P n ~ £(a 5 >/n)*(a, ln(n))' 1 (LL); a,, 1 /V«(NLL); af, a s /y/n, 1 (n (NNLL); . . . 



(1.2) 



for the leading logarithmic (LL), next-to leading logarithmic (NLL) and next-to-next-to leading 
logarithmic (NNLL) order. Prior to the work presented here the RGI bottom mass determination 
from large-« sum rules reached NLL and partly NNLL level [12, 13], but did not include the NNLL 
running of the dominant heavy quark pair production current. Earlier fixed-order analyses [14, 15, 
16, 17, 18, 19] up to next-to-next-to-leading order (NNLO) only resum the Coulomb singularities 
and count (a s \n(n)) 1 as ocj in Eq. (1.2). The convergence of the fixed-order expansion however 
turned out to be rather unsatisfactory, see Ref. [20] for a review. As we will show below RGI 
computations improve the convergence properties substantially and allow for a reliable and precise 
determination of the bottom quark mass from nonrelativistic sum rules. 

For any determination of the bottom mass parameter with a precision at the percent level it is 
mandatory to adopt an appropriate short-distance mass-scheme in order to avoid 0(Aqcd) infrared 
renormalon ambiguities. Suitable mass schemes are the MS scheme for the low-rc sum rules and 
so-called threshold mass-schemes [21, 22, 23, 24] for the large-n sum rules. 

The present talk focuses on the determination of the IS bottom mass [22] from RGI large-« 
sum rules and is mostly based on the recently published Ref. [25]. In this analysis single moment 
fits of the mass parameter are carried out including for the first time the almost complete NNLL 
correction to the theoretical moments. The still missing contribution from the NNLL soft mix- 
ing correction to the running of the heavy quark production current can be neglected under the 
assumption that its size is comparable to the already known soft NNLL terms. 



2. Experimental Moments 

The dominant contribution (87% - 98% for n = 6-12) to the experimental moments in Eq. (1.1) 
for large n comes from the first four T resonances, Y(1S)-Y(4S), which we construct from their 
electromagnetic decay widths and masses [1] using the narrow width approximation. For the con- 
tribution from the threshold region (5.7% - 1.4% for n = 6-12) we use BABAR data [26] in the 
energy range between sfs = 10.62 and y/s = 11.21 and follow the approach of Ref. [27]. Finally 
the continuum region above 11.21 GeV contributes only a tiny fraction. It can be modelled by 
the respective pQCD result [28] assigning a model uncertainty of 10% to the cross section with- 
out introducing a numerically relevant error to the experimental moments. We emphasize that this 
continuum contribution should be regarded as a rough estimate for the (missing) experimental data 
rather than an additional theory input. The precise numbers for the relevant experimental moments 
together with their statistical and systematical errors can be found in Ref. [25]. 
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3. Theoretical Moments at NNLL 

For details on the derivation of the theoretical moments according to the scheme in Eq. (1.2) 
we refer to Refs. [14, 15, 18, 25]. The resummation of nonrelativistic logarithms follows the vN- 
RQCD approach [25]. Here we shall only discuss the general structure and the latest computational 
progress concerning the RGI of the theoretical large-n moments. The theory prediction for the nor- 
malized total bb pair production cross section R^s) in the threshold region is due to the optical 
theorem related to nonrelativistic current-current correlators, which describe the production and an- 
nihilation of a heavy quark pair. Explicit results for these correlators can be adopted from previous 
works on tt threshold production in e + e~ collisions [29, 30, 31, 32]. After the moment integration 
over s the result for the n-th moment through NNLL order can be expressed as 

P " 1 ' NNLL = ~ A l^3S ,n \ Cl & y ) 2 *V & v) + 2c, (h, v)c 2 (h, v) Pn2 (h, v) \ , (3.1) 

where the p„,- arise from the integration of the nonrelativistic current correlators and the c, are 
Wilson coefficients of the respective effective currents. The variables h and v are introduced to 
parametrize the matching and renormalization scales of the effective theory. The natural choice is 
h ~ 1 , V ~ 1 / %fn. The residual dependence of the bottom mass fit on these parameters is used for 
the perturbative error estimate in Sec. 4. 

Equation (3. 1) explicitly depends on the bottom pole mass MP° le , which we translate to the IS 
mass using the relation 

Mf; le = M l b s { 1 + A LL + A NLL + [(A LL ) 2 + a^nll + a nnl L] i _ (3 2) 

The A terms are labeled according to the nonrelativistic order counting scheme in Eq. (1 .2). Explicit 
expressions can be found in Ref. [30]. For the final theoretical expression used in the single moment 
fits below we consistently expand out the perturbative series for the Wilson coefficients c,- together 
with the nonrelativistic expansion series for the p„ ,• and M^ ole in Eq. (3.1). 3 The convergence 
properties of this expansion are discussed in detail in Ref. [25] . 

Apart from the NNLL correction to the renormalization group (RG) running of the Wilson 
coefficient c\ associated with the dominant heavy quark production current all relevant contribu- 
tions to Eq. (3.1) are known completely. Concerning the NNLL running of c\ all ("non-mixing") 
contributions from genuine vNRQCD three-loop diagrams were computed in Ref. [32]. The cor- 
responding pNRQCD calculation is not available at present. The remaining NNLL ("mixing") 
contributions are generated by corrections to the vNRQCD four-quark operator ("potential") co- 
efficients appearing in the NLL anomalous dimension of c\. Recent results for the ultrasoft NLL 
running of the subleading (ff(v) and €?(v 2 )) nonrelativistic quark-antiquark potentials [33, 34, 35] 
completed the calculation of the ultrasoft part of the NNLL mixing contributions [34]. It is the 
dominant NNLL mixing effect [36, 25] in the running of c\. Likewise the ultrasoft terms dominate 
the NNLL non-mixing running [32]. Thus at present the only unknown piece in /^ /! ' NNLL j s the 
NNLL soft mixing contribution to the RG evolution of c\. 

3 We are forced to simultaneously expand out the series for M£ ole and the p„., in the way explained in Ref. [18] in 
order to achieve a proper cancellation of the leading renormalon. 
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Figure 1: RG evolution of the current coefficient c\: NLL (blue) and approximate NNLL result (red) with 
uncertainty due to the unknown NNLL soft mixing contribution (light red band). 



Figure 1 shows the dependence of the current coefficient ci(v) = c\(h = l, v) on the renor- 
malization parameter v for the complete NLL result (blue) and an approximate NNLL result (red), 
where all known NNLL contributions 4 are added. The (light red) band around the NNLL curve is 
generated by varying all known soft NNLL contributions to that curve by a factor between and 
2. The uncertainty due to the unknown NNLL soft mixing terms estimated by this band is much 
smaller than the large total NNLL correction from the running in the relevant range 0.3 < V < 0.6 
and can safely be neglected in the following [25] . 

The (known) leading nonperturbative power correction to the moment p'J 1 ^ 1 ^ j s associated 
with the gluon condensate [14] and turns out to be completely negligible for our analysis [25]. 
Higher order power corrections are sufficiently suppressed for n < 10. 

4. Single Moment Fits 

The vNRQCD expression for the theoretical moment />^< NNLL exhibits a residual dependence 
on the scale /I/, = hnib, where the effective theory is matched to full QCD, as well as on the two 
correlated renormalization scales [1$ = hm^v (soft) and = hmi,v 2 (ultrasoft). Here and in the 
following my = M h ls . The three unphysical scales can be consistently parametrized by the two 
variables h and v. In order to estimate the uncertainties from higher order perturbative corrections 
we choose to vary the parameters for the bottom mass fits within the h-v region around the default 
values v = V* := 1 / y/n + 0.2 and h = 1 as defined in Fig. 2 a. The plot also shows the contours of 
the result for M b ls from the equation P^(M l h s ) = P^f, i.e. a NNLL single moment fit for n = 10. 

With these conventions for the h-v scaling variations we can now generate error bands around 
the default fits (dashed lines) for M b ls (a s ) as shown in Fig. 3. The two panels in this figure compare 
the fits of the bottom mass as a function of the strong coupling a s {Mz) using theoretical moments 
calculated in the fixed-order (a) and the RGI approach (b). The fixed-order moments are obtained 
by switching off the all-order resummation of nonrelativistic logarithms in our RGI moments as 
explained in Ref. [25]. Since at leading order (LO) the only relevant physical scale is the soft scale 

4 These even include the first logarithm °= lnv in the NNLL series of the soft mixing contribution [32]. The 
subset of spin-dependent terms in the NNLL soft mixing contribution is also known [37], but tiny and neglected here. 
All relevant analytic NNLL expressions for c\ are given in Ref. [25]. 
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Figure 2: Panel a): Contour plot of the IS bottom mass determined from P'^my) = P^P as a function of 
the parameters h and / = v/v*. The different contours are labeled by the respective mass value in GeV. The 
region in the h-f plane bounded by the red dashed line represents the parameter space we scan to determine 
the variation of the mass, which contributes to our perturbative error estimate. The region is defined by 
0.75 < h < 1 /0.75 and demanding that 0.5 ju^ < flu < 2 fly, where fi^ = mtv}. The red point inside this 
area indicates our default values f = h=l for the mass determination. Panel b): IS mass results (dots) with 
perturbative error bars from single moment fits for n — 4 to n — 20 as explained in the text. 



the LO and LL bands in Fig. 3 agree exactly. Comparing the next-to-leading (NLO) with the NLL 
and in particular the NNLO with the NNLL results we however observe much larger scale vari- 
ations of the fixed-order results. This clearly indicates a substantially improved precision related 
to the resummation of logarithms in the RGI approach. As argued in Ref. [25] we believe that in 
contrast to the LL band, which is generated only by soft scale variations, and the (w.r.t. the default 
fits) strongly asymmetric NLL error bands the NNLL mass range gives a reliable estimate of the 
perturbative uncertainty. The observed bottom mass dependence on the input value for a s {Mz) is 
rather mild and at least in the interval 0.113 < a s (Mz) < 0.120 linear to a good approximation [25]. 
A far more detailed analysis of the numerical results including plots of Fig. 3 b type for different 
values of n as well as multiple moment fits has been carried out in Ref. [25]. 

The final result for the IS bottom mass from the NNLL RGI single moment analysis for n = 10 
outlined above is 

M l b s = 4.755 ± 0.057 pert ± 0.009a, ± 0.003 exp GeV , (4. 1) 

where we have used the current world average a s {Mz) = 0.1183 ± 0.0010 for the strong cou- 
pling [38]. The central value in Eq. (4.1) corresponds to the result of the default fit (h = 1, V = V*). 
The quoted errors refer to the perturbative uncertainty, which we estimate by half the size of the 
band from the scale variations, and the errors from the uncertainties of a s and the experimental 
data used for the fit. 

Figure 2 b compares the M^ s results and their respective perturbative error bars from fits using 
the moments n = 4 to n = 20. Within the errors all central values are in very good agreement. 
We however emphasize that for n considerably larger than 10 formally the OPE for the theoretical 
moments breaks down due to possibly uncontrolled higher order power corrections, although the 
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Figure 3: Comparison of the masses obtained from the fixed order (a) and RGI calculation (b) of the 10-th 
moment, P&hrib) = P^n '. In panels a and b we show the mass values with LO, NLO, NNLO and LL, NLL, 
NNLL accuracy, respectively. The dashed lines display the results from the fits with the default values for 
the parameters h and v. The corresponding (partly overlapping) error bands were generated by varying h and 
V within the parameter space defined in Fig. 2 a. (We also added the tiny experimental error in quadrature, 
which is however hardly visible.) Concerning panel a, we note that for some low mi, values in the NLO 
band and the associated values for h and v the ultrasoft coupling a s {^iu) reaches 0.65 causing numerical 
instabilities. 

leading power correction still appears to be small as long as n < 20. We therefore regard the 
error bars for n > 10 shown in Fig. 2 b as a confirmation of the perturbative stability of the RGI 
vNRQCD calculation, but do not use them for quoting our final errors. On the other hand the error 
bars for n < 10 increase for smaller n because the sensitivity of the theoretical moments on the 
mass decreases [25], cf. Eq. (3.1). 

Using the respective (fixed-order) relations to the pole mass, see Refs. [18, 19] for details, we 
can translate the IS mass result in Eq. (4. 1) to the MS-mass and obtain 

m b {m b ) = 4.235 ± 0.055pe rt ± 0.003 exp GeV , (4.2) 

where we have added an additional conversion error of 15 MeV to the perturbative uncertainty [25]. 
Interestingly the a s dependence of the original IS mass result in Eq. (4.1) and the intrinsic a s 
dependence of the IS -MS conversion formula almost cancel exactly. The remaining a s induced 
error is therefore negligible and not quoted in Eq. (4.2). 

5. Summary 

We have presented the determination of the IS bottom mass from (single) large-« moment 
fits with RGI at NNLL order as carried out in Ref. [25]. The main result is given in Eq. (4.1). 
Converted to the MS scheme our result (Eq. (4.2)) is consistent with the NLL RGI large-« result of 
Ref. [12], but not quite compatible with the latest results from low-« sum rule determinations [3, 5]. 
We however note that our calculation (like Ref. [12]) treats the charm quark as massless, while 
previous fixed-order analyses [39, 19] have shown that finite charm mass effects are enhanced for 
large n and cause a sizable mass shift between —20 and —30 MeV. A similar effect is expected in 
the RGI analysis and might help to reconcile the discrepancy. 
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